The protein dynamical transition is a pseudogap changeover 
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The emergence of biochemical activities in a protein seem to commence with the onset of atomic 
mean-square displacements along the protein lattice. The ensuing protein dynamical transition has 
been discussed extensively in the literature, and often with conflicting conclusions. Here we clarify 
the phenomenon by establishing a deep connection between the dynamical transition and the pseu- 
dogap state where high-temperature superconductivity comes to its end. For this we first show how 
to endow proteins with an order parameter akin the quasiparticle wave function in superconductors. 
We then present universality arguments to claim that the protein dynamical transition takes place 
in tandem with a pseudogap transmutation. We confirm that available experimental data fully 
supports our proposal. 



The protein dynamical transition was first observed 
around 30 years ago in myoglobin ^-[3]. It is now 
claimed to be a common property of all hydrated proteins 
[4]- [6]. The transition takes place at temperatures that 
are somewhere between 180K-240K. Since this coincides 
with the temperature range where proteins quite univer- 
sally start to display measurable biochemical activities, 
it has been proposed that the protein dynamical transi- 
tion is concomitant to the onset of life |4]-[6]. But the 
detailed nature of the protein dynamical transition, even 
its very existence, remains controversial [7] , [H] . The phe- 
nomenon, if it indeed exists, appears to closely mimic the 
properties of the supercooled water that surrounds the 
protein [9], [10]: At temperatures below ^ 150K a flash- 
cooled water under ambient pressure is in an amorphous, 
glassy state. Between ~ 150K-240K there is a no man's 
land where the time scale for crystallization is too short 
to be analyzed with present day techniques. But there is 
evidence of an ice-liquid coexistence in this regime |10j . 
At around 240K water then enters a crystalline state with 
a homogeneous ice nucleation that persists until the liq- 
uid phase takes over at 273K. 

To a large extent protein folding is known to be driven 
by a combination of hydrophobic and hydrophilic effects. 
Therefore the properties of a hydrated protein should 
mirror those of the adjacent water. Indeed, there is a 
wide consensus that the protein dynamical transition is, 
likewise a glass transition, driven by two intertwining 
processes that are controlled by a combination of the dif- 
ferent properties that water has in the thin hydration 
shell and in the surrounding bulk. At an atomary level 
the two processes are as follows |4]-[6]: At temperatures 
below ^ 180K the covalent protein lattice is in a state 
where the individual atoms are predominantly subject to 
local thermal fluctuations. These fluctuations can be in- 
terpreted as simple harmonic vibrations of atoms around 
their lattice equilibrium positions. The Debye- Waller re- 
lation connects the atomic mean square displacements 
to the experimental B-factors, that grow linearly at low 
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When the protein dynamical transition takes place, 
atomic displacements begin to correlate and start cov- 
ering much larger length scales. Various collective mo- 
tions such as fluctuations between different macromolec- 
ular sub-states make the scene, take over and rule the in- 
ternal dynamics of the protein lattice all the way to phys- 
iological temperatures. This becomes reflected in the B- 
factors that should now display a more rapid, maybe even 
exponential increase as a function of temperature. 

The conventional approach to the protein dynamical 
transition is very much concentrated on understanding 
the fine structure and detailed properties of supercooled 
water 4j-[6 . Here we propose a radical departure from 
the paradigm way. Our approach is based on the theoret- 
ical concept of universality that already has a firm basis 
in the description of phase properties of matter in terms 
of a relatively small number of relevant and marginal in- 
teractions [llj . We employ a combination of universality 
arguments with experimental data analysis to conclude 
that the protein dynamical transition must be a pseudo- 
gap changeover. We argue that the pertinent pseudogap 
state is the one that precedes the 0-point phase transi- 
tion where proteins become denatured and depart from 
their biologically active collapsed phase. Furthermore, 
since our arguments are entirely based on the concept of 
universality, the transition can not be specific to proteins 
and water only, it is a more general property of polymers 
in bad solvents. 

A pseudogap state was originally introduced to explain 
aspects of high temperature superconductors [T^], [13] ■ 
Subsequently it has been found experimentally even in 
ferromagnetic metals [H]. Theoretically, a pseudogap 
has been proposed to be a participant both in color su- 
perconductors and in the chiral phase transition of QCD 
[Is] . [16] . At a theoretical level and in its simplest form, 
the pseudogap state is described by a single complex or- 
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der parameter (f) = p exp{i0}. In the case of superconduc- 
tivity this could be the order parameter for the Cooper 
pairing of electrons. In the symmetric phase where there 
is no superconductivity, the free energy has its minimum 
when the modulus < >=< p > vanishes. The symme- 
try becomes broken when there is condensation of (f> and 
the free energy is minimized by a non- vanishing value of 
the modulus < p >. In the case of a superconductor this 
takes place when electrons combine into Cooper pairs, 
the hallmark of a superconducting phase. The pseudo- 
gap state is a refinement of this phase structure into the 
case where < 4> > vanishes even though < p > retains its 
non- vanishing symmetry breaking value |12j . This occurs 
in the presence of a strong phase decoherence. 
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Consequently we can have < >= even though the 
nonvanishing gap < p > persists. The pseudogap state 
is like a symmetric phase precursor state in the broken 
symmetry phase, a prelude to the transition that takes 
place when the gap < p > eventually vanishes. 

In the case of proteins, more generally of polymers, the 
conventional order parameter is the compactness index v 
that describes how the radius of gyration Rg scales in the 
limit where the number N of residues becomes very large 
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Here (i ~ 1,2,...,N) are the locations of the central 
Cq carbons and Rq is a protein specific but N indepen- 
dent pre-factor. Canonically, the compactness index can 
assume only three different values. The mean-field the- 
ory value v = 1/3 corresponds to the biologically active 
collapsed phase where proteins are in a space filling con- 
formation. This commonly occurs at low temperatures 
or in a bad solvent environment. When v = 1/2 we have 
a fully flexible chain, this corresponds to the 0-point. 
Finally, = 3/5 is the mean-field Flory value for the 
universality class of self-avoiding random walk that de- 
scribes a polymer either at a very high temperature or 
in a very good solvent environment. This conventional 
approach to protein phase structure does not have the 
flexibility to describe a pseudogap state, so we now pro- 
pose a reflnement. We start by introducing the local 
backbone bond and torsion angles. We utilize the Cq 
coordinates to deflne the backbone unit tangent (t) and 
binormal (b) vectors 
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Then 



are the discrete bond and torsion angles, respectively. 
Note that these deflnitions involve only the positions of 
the Cq carbons. Consequently each quantity is inherently 
geometrical and we can introduce the similarly geomet- 
rical order parameter 
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In fact, this is exactly the Hasimoto identification of the 
wave function in nonlinear Schodinger equation |18| . in 
its discrete version. As a complex variable the order pa- 
rameter ([6]) also has the requisite structure to describe a 
pseudogap. We now proceed to argue that in a collapsed 
protein the pseudogap state also occurs, and that it is a 
precursor state to the 8-point phase transition. 

In [19J it has been shown that most proteins in Pro- 
tein Data Bank, over 90% of them, can be described with 
experimental B-factor accuracy in terms of soliton solu- 
tions to the following variant of the discrete nonlinear 
Schrodinger equation [20] 
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The first sum together with the three first terms in the 
second sum comprise exactly the energy of the standard 
discrete nonlinear Schrodinger equation (DNLS) when 
expressed in terms of the Hasimoto variable ([6]). The 
fourth {br) and fifth [a-r) terms are the two conserved 
quantities that precede the energy in the DNLS hierar- 
chy, the momentum and the helicity |21j. The last (c^.) 
term is the Proca mass term, we include it only for com- 
pleteness. 

The energy function ([t]) does not purport to explain 
the details of the atomary level mechanisms that give 
rise to protein folding. Rather, it allows us to exam- 
ine the properties of a folded protein backbone in terms 
of universal physical arguments, much like an effective 
Landau- Ginzburg model describes superconductivity. In 
fact, we can further develop the analogy between ([6]) and 
Cooper pairs, by noting that ([?]) has the functional form 
of the discretized Landau-Ginzburg free energy in the su- 
percurrent variables |i22j : In the continuum limit the first 
two terms combine into the derivative of k{x) that plays 
the role of Cooper pair density. The third term is the 
standard symmetry breaking potential and the fourth 
term has its origin in spontaneous symmetry breaking 
that leads to the notorious Meissner effect [12] ■ In addi- 
tion we have included in ([t]) exactly all those terms that 
are consistent with general principles of universality and 
gauge invariance. 

If we use the Ti equation of motion to eliminate the 
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torsion angles, the potential for the bond angles becomes 
^ ad + be „ / 9 b'^ \ 9 4 /„N 

The first term relates to the potential energy in a 
Calogero-Moser system, notorious for its role in describ- 
ing fractional statistics. Here it is slightly deformed by 
the Proca mass. The second and third terms have the 
conventional form of a symmetry breaking double-well 
potential. Depending on the parameter values we may 
be either in the symmetric Ki = phase or in the bro- 
ken symmetry phase where Ki acquires a non-vanishing 
ground state value. In the second case the local en- 
ergy minimum states correspond to regular protein sec- 
ondary structures such as a-helix or /3-strand and pro- 
tein loops are like domain wall configurations that inter- 
polate between the different local minima [53]. Explicit 
computations confirm |24j that the symmetry breaking 

« 1/3 low temperature phase where long-range corre- 
lations rule, becomes converted into a symmetry restored 
ly ~ 1/2 phase at higher temperatures. In particular, 
the symmetry restoration is a phase transition that takes 
place when temperature reaches the 0-point value that 
separates the low temperature collapsed phase from the 
higher temperature ideal chain phase [24] . 

Unlike in the case of conventional metallic supercon- 
ductors where the supercurrent has stiffness that leaves 
no room for a pseudogap, the torsion angles in ([t]) are 
only subject to ultra-local couplings and as such they are 
highly flexible. Thus it is conceivable that as temperature 
starts increasing in the low temperature phase, there will 
also be an increase in phase decoherence that is measured 
by < exp{ir} >. Depending on the dynamical details, 
the energy function ([t]) may then describe a pseudogap 
state at temperatures below the 0-point value. But as it 
stands, the present energy function can not fully model 
the pseudogap transition in proteins. This is because the 
model does not include the side-chains with their numer- 
ous rotamers Xi- The rotamers give us plenty of opportu- 
nities to introduce many additional pseudogap detecting 
order parameters < expjzxi} > in proteins: 

At low temperatures protein is in a crystalline state 
and the rotamers assume very definite values that are cat- 
alogued in various libraries. As temperature raises fluc- 
tuations appear, but these are damped by various steric 
constraints that strongly limit the ability of individual 
atoms to move. In proteins the steric constraints are 
more lenient to the outlying side-chain rotamers than to 
the backbone torsion angles: With fever covalent neigh- 
bors, there is more room to fluctuate. Consequently any 
protein pseudogap transition should become primarily 
visible in an increased phase de-coherence in the outlying 
side-chain rotamers. But even though steric constraints 
are more lenient to side-chain rotamers, the identification 
of a protein pseudogap is much more delicate than the 
identification of a pseudogap state in materials such as 



high-temperature superconductors where no steric con- 
straints exist. Despite being numerous, the protein side- 
chain rotamers can never assume unhindered, mutually 
fully uncorrelated values. Consequently protein rotamer 
angles can never display a full, unrestrained phase de- 
coherence the way how it is laid out by a high temper- 
ature superconductor. Rather we expect there to be a 
pseudogap changeover that makes its presence known in 
a marked, rapid relative increase in the side-chain ro- 
tamer fluctuations as the temperature climbs up towards 
the 9-point value. 

Presently, there are no direct experimental measure- 
ments of temperature dependence in the rotamer pseu- 
dogap order parameters < exp{ixi} >■ We propose that 
such experiments could be performed. At the moment 
our ability to investigate the inception of a pseudogap is 
limited to the analysis of increased temperature depen- 
dence in the experimental side-chain B-factors that indi- 
rectly measure the < expjjXi} > fluctuations. We have 
screened all data available in Protein Data Bank, but mu- 
tually compatible good quality PDB data on B-factors is 
sparse. Structures have been determined without stan- 
dardized crystalline environments, using different refine- 
ment procedures, and with apparently different practices 
for determining the B-factors. This makes it very hard 
for us to compare data that is collected in different ex- 
periments and at different temperatures. Consequently 
our ability to reliably investigate the relation between a 
pseudogap changeover and the protein dynamical tran- 
sition is limited to the re-analysis of data that has been 
collected in only two experiments, with crambin [25] and 
ribonuclease A [5S] . They both measure B-factors at dif- 
ferent temperatures, consistently with the same protein 
crystals. We have re-analyzed the data and find that both 
experiments fully support our proposal that the protein 
dynamical transition is a pseudogap changeover. 

In Figure 1 we display the B-factors of these two exper- 
iments as a function of temperature, separately for the 
backbone C^ carbons and for the side-chain atoms. From 
the side-chains we exclude the C^ atoms, they are slaved 
to the sp3 hybridized backbone Cq atoms both by cova- 
lent bonds and by existing refinement procedures. In the 
case of crambin we observe a clearly visible qualitative 
change and a definite enhancement in the temperature 
dependence of the side-chain B-factors in comparison to 
the backbone Cq atoms. There is also a transition from a 
linear to a non-linear behavior. The change commences 
with the protein dynamical transition regime. For ri- 
bonuclease A the results are similar, but the effect is less 
profound. We remark that in both cases the B-factors of 
the Cq atoms truly deviate only slightly (if at all in the 
case of ribonuclease A) from linearity in T. This is con- 
sistent with [27] . where no sign of a protein dynamical 
transition was observed at the level of the backbone Cq. 
carbons. 

In summary, the protein dynamical transition and the 
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FIG. 1: a) The B-factors of crambin from [25| as a function of 
temperature, separately for backbone (red) and side-chain 
atoms (black) but excluding the C^. b) Same as a) but for 
ribonuclease A, with data from [2^. In both cases the side- 
chains display a clear transition and nonlinear T dependence, 
while for the Cq even the linear fit ([TJ is feasible. 



pseudogap state in high temperature superconductors are 
two a priori totally different physical phenomena. How- 
ever, here we have proposed that they are intimately re- 
lated by asserting that the protein dynamical transition 
is a pseudogap changeover. Our argumentation is based 
both on theoretical observations that stem from the con- 
cept of universality, and analysis of all presently available 
experimental data. We propose that experiments could 
be performed to directly detect the pseudogap transition, 
in particular since the possibility that a pseudogap con- 
verts a static protein crystal into a biologically active 
nano-engine could have really far reaching consequences 
in the search for a physical origin of life. 
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